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ABSTRACT 

We present the orbital and physical parameters of a newly discovered low-mass 
detached eclipsing binary from the All-Sky Automated Survey (ASAS) database: 
ASAS J011328-3821.1 A - a member of a visual binary system with the secondary 
component separated by about 1.4 seconds of arc. The radial velocities were calcu- 
lated from the high-resolution spectra obtained with the 1.9-m Radcliffc/ GIRAFFE, 
3.9-m AAT/UCLES and 3.0-m Shane/HamSpec telescopes/spectrographs on the ba- 
sis of the TODCOR technique and positions of H a emission lines. For the analysis we 
used V and / band photometry obtained with the 1.0-m Elizabeth and robotic 0.41-m 
PROMPT telescopes, supplemented with the publicly available ASAS light curve of 
the system. 

We found that ASAS J011328-3821.1 A is composed of two late- type dwarfs having 
masses of Mi = 0.612 ± 0.030 M , M 2 = 0.445 ± 0.019 M Q and radii of i?i = 0.596 ± 
0.020 R Q , i? 2 = 0.445 ± 0.024 R Q , both show a substantial level of activity, which 
manifests in strong H a and Hp emission and the presence of cool spots. The influence of 
the third light on the eclipsing pair properties was also evaluated and the photometric 
properties of the component B were derived. Comparison with several popular stellar 
evolution models shows that the system is on its main sequence evolution stage and 
probably is more metal rich than the Sun. We also found several clues which suggest 
that the component B itself is a binary composed of two nearly identical ~ 0.5 M Q 
stars. 

Key words: binaries: eclipsing - binaries: spectroscopic - binaries: visual - stars: 
fundamental parameters - stars: individual (ASAS J010538-3821.1) - stars: low-mass 



1 INTRODUCTION 

In the last decade there has been a stunning increase in 
our knowledge of the structure and evolution of eclipsing 
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binaries. The number of known and well studied detached 
low-mass eclipsing systems, which are the main source of 
absolute parameter measurements, increased by almost an 
order of magnitude after the year 2000. This increase helped 
several authors make conclusions about a 30-year old prob- 
lem: the discrepancy between observed and theoretically pre- 
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dieted radii and temp eratures during the evolution of low- 
mass stars l|Lacvll 19771 ). The prevailing idea is that this dis- 
crepancy is caused by the stellar magnetic field, amplified 
by the fast rotation due to tidal locking in close binaries, 
which inhibits the efficiency of convection in the envelope 
and manifes ts itself in a higher level of activity than for 
single stars l|Chabrier et al.|[2007l ). However, the number of 
known well-characterized systems is too low to reliably test 
this hypothesis. This is especially true amoung the long- 
period systems, which normally exhibit lower levels of ac- 
tivity, so the observed radii and temperatures are closer to 
the theoret i cally p redicted values. In a very recent paper, 
llrwin et alj l|201ll ) described an eclipsing pair of M dwarfs 
on a 41 d orbit which, despite slow asynchronous rotation 
and a probable age over 120 Myr, still show cold spots, and 
radii are inflated by 5% with respect to theoretical predic- 
tions. This system is evidence against the theory mentioned 
above and may again bring wider interest to the need to 
revisit the equation of state for low mass stars, which was 
suggested to ex plain the inflated rad ii before the activity 
hypothesis (e.g. (Torres fc Ribadl200l ). 

To be useful for testing the evolutionary models, the pa- 
rameters of the system, especially masses and radii, should 
be known with a pr ecision of at least 3%, considered now 
as a canonical level (jBlake et al. 2008; Clausen et al.l l2008). 
There are two other criteria which can make a system more 
interesting: (1) additional multiplicity - by deriving prop- 
erties of the third body one can put additional constraints 
on the nature of the whole system; (2) low mass ratio - in 
general it is harder to fit a single isochrone to data points 
that are well separated in parameter space; these are usually 
stars of different masses. Twin stars commonly have data 
points that lie close together and are less difficult to fit. 
In this paper we present our analysis together with orbital 
and physical parameters of a newly discovered low-mass de- 
tached eclipsing binary which meets the two mentioned cri- 
teria. The data we gathered fell short of the 3% precision 
level goal, but the system is still highly valuable for our 
understanding of the nature of low-mass stars. 



2 THE ASAS J011328-3821.1 SYSTEM 

The eclipsing system ASAS J011328-3821.1 (2MASS 
J01132817-3821024, 1RXS J011328.8-382059; hereafter 
ASAS-01) has the shortest orbital period from all ob- 
jects in our sa mple of low-mass detached eclipsing binaries 
(LMD EBs; see: lHelminiak fc Konackill201ll ; lHelminiak et all 
2011) f ound in the AS AS Catalogue of Variable Stars 
(ACVS; |Poimanskii r2002fl. and one of the shortest among 
the LMDEBs known to date - Pasas = 0.44559 d. In the 
ACVS the ASAS-01 is classified as eclipsing detached (ED) 
with a maximum brightness of V = 11.78 mag and ampli- 
tude of brightness variation of 0.33 mag. The ASAS light 
curve clearly shows that there is a significant difference in 
the eclipse depths which implies that the eclipsing pair is 
composed of two vastly different stars. 

Before the ACVS the object was not reported to be 
a binary or multiple but was noted as an X-ray source 
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Figure 1. A part of an image taken with 1.0-m Elizabeth tele- 
scope in SAAO showing the two components of the ASAS-01 AB 
system. The visible components are separated by ~1.25 asec. 
ASAS-01 Aa + Ab is a low-mass eclipsing detached binary identi- 
fied in the ASAS database (Pas AS = 044559 d). The orientation 
and the scale of the image are shown. 



in the Ei nstein dGioia et al.l Il990l ; IStocke et alj 1991 ) and 
ROSA T (|Voges et all 1 19991 ) catalogues. Later iBeers et all 
jl996l ) noted emission in calcium H and K lines, which they 
describe d in their c a talogu e as "moderate to strong". More 
recently iRiaz et al.l (120061 ) identified ASAS-01 as an M0. 5 
dwarf, gave values of spectral indexes CaHl, CaH2, CaH3 
and Ti05, reported H a emission line with equivalent width 
of 2.5 A, and estimated the distance to be about 37 pc. 
This determination is significantly diff e rent f rom the value 
of 26.6 pc obtained by ISzczygiel et af] (|2008t) . who derived 
it on the basis of a bolometric correction derived from the 
ASAS V — I colour, but the object was treated as a single 
star, explaining the discrepancy. 

The emi s sion in H a was recently confirmed by 
iParihar et alj |2009) who estimated the spectral types of 
components to be M1V + M3V. Their study, however, did 
not reflect the possibility of the system being composed of 
3 or more stars. 

Finally, in their lucky-imag ing survey of southern 
M dwarfs, iBergfors et alj f2010) obtained high-angular- 
resolution images of the system in i' and z bands. They 
clearly separated two components, with the secondary lo- 
cated 1.405(3) seconds of arc at the position angle of 29.0(3) 
degrees. They estimated the spectral types to be M0 + Ml, 
not taking into account the eclipsing nature of one of the 
stars. 

We deduced the multiple nature of ASAS-01 in two 
ways. Firstly we noted that in the spectrum taken with the 
Shane telescope in 2007 (Sect. the H a emission line profile 
was found to be triple. The same feature was found in a spec- 
trum obtained in 2008 with the Anglo- Australian Telescope 
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but not in data obtained in 2006 with the Radcliffe tele- 
scope in South Africa . Shortly thereafter and independent of 
iBergfors et ail l|2010h we found the two components clearly 
separated on the CCD images obtained with the Elizabeth 
telescope in South Africa (Sect. |4} in good seeing conditions. 
A part of our image is shown in Fig. [T] The fainter B compo- 
nent is located about 1.25 arcseconds from the brighter A in 
the NE direction, which is consistent w ith the more precise 
result obtained by IBergfors et ail (|2010l ). A series of images 
taken in different orbital phases (with respect to the photo- 
metric 0.44559 d period) allowed us to deduce that eclipses 
occur in the component A. 



3 SPECTROSCOPY AND RADIAL 
VELOCITIES 

3.1 Observations 

To directly derive masses of ASAS-01 components we ob- 
tained a series of high-resolution spectra in order to mea- 
sure their radial velocities (RV). Most of the spectroscopic 
observations come from the 1.9-m Radcliffe telescope and its 
Grating Instrument for Radiation Analysis with a Fiber-Fed 
Echelle (GIRAFFE) spectrograph at the South African As- 
tronomical Observatory (SAAO), with sett ings the same as 
descr i bed i n previous papers of this series (|HeIminiak et all 
120091 1201 ll ). Five spectra of R ~ 40000, were obtained dur- 
ing three consecutive nights in September 2006. After 3600 
s of exposure and binning to R ~ 10000 we obtained SNR 
at A = 6520 A between 25 and 40. At that time we were 
not aware that the system is also a visual binary, and the 
quality of the image from the aquisition camera was not 
enough to demonstrate that. Due to a relatively high seeing 
(about 2 asec and more) our spectra were contaminated by 
the companion B. 

One additional spectrum was obtained in October 2007 
with the 3.0-m Shane telescope and its Hamilton Spectro- 
graph (HamSpec) at the Lick Observatory. The target was 
observed very low over the horizon for 2100 seconds. After 
binning from the original resolution of 60000 to 30000 we 
reached the SNR level around 25 at A = 6560 A. Due to 
the low SNR, fringing, telluric lines and limited extent of 
available templates we did not use all 96 rows of the spec- 
trum, but only used 21 that covered the wavelength range 
from 5150 to 6500 A. This spectrum was also contaminated 
by the light from the component B. 

The last spectrum was obtained in September 2008 
with the 3.9-m Anglo- Australian Telescope and its Univer- 
sity College London Echelle Spectrograph (UCLES) at the 
Siding Sp ring Observatory , with the s ame settings as de- 
scribed in Hclminia k et al.l l|2009l . l2~01lf ). This spectrum was 
not binned and after 900 s of exposure we reached SNR ~ 20 
at A = 6530 A. We tried to keep the slit on component A 
only, by rotating it by 90 deg from the position angle of the 
AB pair, but due to the high seeing (over 2 asec) we could 
not avoid a significant contamination from star B. 



3.2 Radial velocities derivation. 

The CCD reduction, spectrum extraction and wavelength 
calibration were done in every case with standard iraf pro- 



cedurefl For the wavelength calibration we used exposures 
of ThAr lamps taken before and after the science exposure. 
Initially the RV's were measured with our implementation 
of the two dimensional cross-correlation technique (todcor; 
IZucker fc MazeblH994h using a spectrum of an M2V stan- 
dard star HIP 93101 broadened to w rot sini = 30 km s _1 
that was taken as a template. The formal RV measurement 
errors were computed from the bootstrap analysis of TOD- 
COR maps created by adding randomly selected single-order 
TODCOR maps. 

In the TODCOR maps, the location of the maximum of 
the cross-correlation function in the t>i/v2 plane refers to ra- 
dial velocities of the two components of the system. We could 
easily phase-fold our measurements of the primary t>i with 
the photometric period of 0.44559 d, but we failed when try- 
ing to do this for V2 (putative secondary). Thus to obtain the 
RV curve of the eclipsing secondary we decided to measure 
positions of the H a emission lines. In all GIRAFFE spectra 
we could easily recognize two well-separated components of 
that line and measure their positions despite that the line 
was close to the edge of an echelle order. To measure the po- 
sitions we used the deblend procedure in the iraf/echelle 
package and fitted a double- gaussian function. The fit was 
performed about 20 times with different initial wavelength 
range selections, initial peak positions and smoothing fac- 
tors. Results were then averaged and the standard deviation 
was taken as the measurement error. A similar procedure 
was applied to the HamSpec and UCLES spectra but with 
a triple-gaussian function, since the triple character of the 
H a line was obvious. We did not normalize the spectra to 
the continuum before this procedure. 

To check whether the RV's obtained from H a lines can 
be trustful we tried to phase-fold them separately with the 
photometric period. For the primary we found a good agree- 
ment with the TODCOR measurements, and for the secondary 
we succeeded to find a fit with the desired period, systemic 
velocity close to the one found for the primary, and ampli- 
tude Ki significantly higher than the primary's K\. This 
is expected for two distinctive late-type dwarfs, consider- 
ing the difference in the depths of the eclipses in the ASAS 
light curve. For the final orbital fit we used all measurements 
simultaneously. We used the spectroscopic orbit fitting pro- 
cedure described in our previous papers ([H clminia k et al.l 



2009; iHelminiak fc Konackil l201ll ; iHelminiak et all 12011 



This simple code uses a Levenberg-Marquardt minimaliza- 
tion algorithm to find a keplerian orbit of a spectroscopic 
binary. It also allows for estimation of systematic contri- 
butions to the error budget by Monte-Carlo and bootstrap 
analysis. 



3.3 ASAS-01 A orbital fit 

Before obtaining the orbital fit we used our photometric 
measurements (Sect. [4} to correct the ASAS ephemeris. Us- 



2 IRAF is written and supported by the IRAF programming group 
at the National Optical Astronomy Observatories (NOAO) in 
Tucson, AZ. NOAO is operated by the Association of Uni- 
versities for Research in Astronomy (AURA), Inc. under co- 
operative agreement with the National Science Foundation, 
http : //iraf .noao . edu/ 
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Figure 2. Left: Radial velocity measurements and the best-fitting orbital solution for the ASAS-01 A eclipsing binary. Open symbols 
and the solid line are for the primary, filled symbols and the dashed line are for the secondary. Circles denote TODCOR , stars denote 
H a lines measurements and crosses mark measurements rejected from the final fit. Residua with the final error bars are shown on the 
lower panel; Right: Radial velocity measurements of the ASAS-01 B component from the TODCOR technique and the best fit of a putative 
circular orbit of P ~ 0.475 d period. Crosses mark measurements rejected from the final fit. Residua and error bars are shown on the 
lower panel. 



ing the jktebop and phoebe codes (see Sect. [S]) we found 
a new value of the orbital period P = 0.44559604(18) d and 
incorporated it to the orbital fit. We have also found (from 
the photometry and radial velocities) that the eccentricity 
e agrees with zero well within formal errors, so we held it 
fixed later. We found the same for the difference between 
systemic velocities of the primary and secondary 7Aa — 7Ab, 
and for any possible systematic shifts between vi measure- 
ments from TODCOR and H a lines. However, the reduced \ 2 
of the best-fitting solutions made separately for H a mea- 
surements of either the primary or the secondary were not 
close to 1. This means that the measurement errors were un- 
derestimated, so we added in quadrature a systematic term 
of 6.5 km s -1 to the formally derived measurement errors. 
This operation was not necessary for the TODCOR results. We 
also excluded one H a measurement of the primary from the 
GIRAFFE spectrum at JD~2454008.477 and both from the 
UCLES spectrum, due to their very large formal uncertain- 
ties. The measurements from the UCLES spectrum probably 
also suffered from the line blending. For this exposure the 
velocity difference was the smallest in the whole sample and 
the emission from the B component is strongest. However, 
we believe that the TODCOR measurement is secure because 
the peak was very high and well separated from other smaller 
maxima. 

All the RV mesurements for all components of the 
ASAS-01 multiple system are listed in Tabled] For ASAS- 
01 A we show only the measurements which were taken to 
produce the orbital fit and the final model of the system, 
together with their final errors and residuals. For the com- 
ponent B we show all the TODCOR V2 values and two H a ve- 
locities derived from HamSpec and UCLES spectra. In the 
"T/S" column "R/G" denotes Radcliffe/GIRAFFE, "S/H" 
denotes Shane/HamSpec and "A/U" denotes AAT/UCLES 
telecope/spectrograph. 



In the left panel of Fig. [5] we present all the measure- 
ments together with the best fitting orbital solution obtained 
for the ASAS-01 A eclipsing pair. We found the solution with 
a fast procedure which fits a double-keplerian orbit by mini- 
mizing the x 2 function with a Levenberg-Marquard method. 
The rms of the fit was 5.8 km s _1 for the primary and 7.0 
km s _1 for the secondary. The final reduced \ 2 was 0.988, 
thus we believe that the parameter uncertainties are well 
estimated. To ensure the accuracy of the estimates we also 
performed an additional Monte-Carlo simulation to estimate 
the systematic errors of derived parameters. The resulting 
systematic terms were an order of magnitude smaller than 
the formal ones, but were added in quadrature to the errors 
from the least-squares fit. 

The orbital parameters derived from the fit are listed in 
Table [5] One can see that H a and TODCOR measurements 
are complementary and, when combined, allow us to reach 
about 2% level of precision in the velocity amplitudes if 1,2 
and about 4% in Msin 3 i. We would like to emphasize the 
low value of the mass ratio q = 0.727(18), which makes 
ASAS-01 A quite unique. With this q ASAS-01 A resides 
in the mass-ratio distribution of spectroscopic binaries well 
outside both th e narrow pe ak at q > 0.95 (the "strong" 
twin hypothesis; lLucvll2006t ) and the wide one at q > 0.85 
(the "weak" twin hypothesis; lHalbwachs et al.l 12003). This 
makes ASAS-01 A particularly interesting for testing the 
evolutionary models. 

3.4 The component B 

The positions of the third H a emission peak in the Ham- 
Spec and UCLES spectra were close to (however not in good 
agreement with) the corresponding measurements of V2 from 
TODCOR . This led to a conclusion that the V2 coordinate of 
the TODCOR maxima may be related to the B component, 
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Table 1. Single RV measurements of all components of the ASAS- 
01 system derived with the TODCOR technique and from the po- 
sition of H a emission lines, together with their formal errors and 
final orbital fit residua. "R/G" denotes data points from Rad- 
cliffc/GIRAFFE observations, "S/H" from the Shane/HamSpec 
spectrum and "A/U" from AAT/UCLES. For the B component 
the residua refer to the fit of a circular orbit with Pg ~ 0.475 d. 
When O — C is not given, the measurement was not included into 
the fit. 



JD 


V 


± 


O 


- C 


T/S 


-2450000 


(km s ) (km s 1 ) (k: 


III £ 


- 1 ) 




Component 


Aa velocities from TODCOR 








4006.50332 


-86.882 


4.360 


3 


.502 


R/G 


4007.52761 


98.787 


7.240 


3 


.361 


R/G 


4007.57128 


140.808 


7.037 


7 


.444 


R/G 


4008.47697 


140.128 


4.695 


3 


.230 


R/G 


4008.52042 


112.337 


4.540 


-6 


.065 


R/G 


4375.87729 


-100.752 


8.175 


-1 


.146 


S/H 


4727.14578 


68.250 


1.247 





.076 


A/U 


Component 


Aa velocities from H a lines 








4006.50332 


-90.330 


6.971 





.054 


R/G 


4007.52761 


96.905 


7.435 


1 


.479 


R/G 


4007.57128 


132.301 


6.895 


-1 


.063 


R/G 


4008.47697 


122.606 


7.232 




— 


R/G 


4008.52042 


107.616 


6.577 


•10 


.786 


R/G 


4375.87729 


-110.900 


6.563 


-11 


.295 


S/H 


4727.14578 


94.415 


6.540 






A/U 


Component 


Ab velocities from H a lines 








4006.50332 


161.590 


7.509 


-7 


.093 


R/G 


4007.52761 


-85.570 


7.401 


1 


.556 


R/G 


4007.57128 


-144.540 


7.134 


-5 


.276 


R/G 


4008.47271 


-136.700 


7.866 


6 


.691 


R/G 


4008.52042 


-113.824 


7.176 


1 


.881 


R/G 


4375.87729 


191.210 


6.655 


9 


.812 


S/H 


4727.14578 


-68.139 


6.525 






A/U 


Component 


B velocities from 


TODCOR 








4006.50332 


116.254 


6.240 


3 


.502 


R/G 


4007.52761 


28.122 


1.500 


3 


.361 


R/G 


4007.57128 


-21.130 


6.076 


7 


.444 


R/G 


4008.47697 


19.200 


3.690 


3 


.230 


R/G 


4008.52042 


-24.688 


9.086 


-6 


.065 


R/G 


4375.87729 


-12.785 


5.925 


-1 


.147 


S/H 


4727.14578 


21.381 


0.482 





.076 


A/U 


Component 


B velocities from 


H a lines 








4375.87729 


13.147 


6.587 






S/H 


4727.14578 


6.799 


6.633 






A/U 



Table 2. Orbital parameters of the ASAS-01 A eclipsing binary. 



Parameter 


Value 


± 


Pa [d] 


0.44559604 


1.8e-7 


T 0jA [JD-2450000] 


1868.83933 


3.1e-4 


K Aa [km s" 1 ] 


118.4 


2.0 


K Ab [km s _1 ] 


162.9 


3.3 


[km s _1 ] 


18.6 


2.3 


1A 


0.727 


0.018 


a a sini [Rq] 


2.478 


0.034 


e A 


0.0 


(fix) 


M Aa sin 3 i [M ] 


0.595 


0.027 


M Ab sin 3 i [M©] 


0.432 


0.017 



Note: Values of P, Tq arc taken from the simultaneous spectroscopic and 
photometric analysis performed with PHOEBE. 



Table 3. Parameters of the ASAS-01 B putative orbit. 



Parameter 


Value 


± 


Pb [d] 


0.4750 


0.0009 


T , s [JD-2450000] 


4806.950° 


0.022 


K Ba [km s _1 ] 


104.7 


5.7 


7 S [km s" 1 ] 


13.6 


2.9 


es 


0.0 


(fix) 


a,Ba sini [R ] 


0.983 


0.053 


/(Affl«) [Mq] 


0.0565 


0.0092 



a Tq is calculated for the moment of the orbital conjunction. 

not the secondary component of the eclipsing Aa+Ab pair. 
If so, the B component itself must be a spectroscopic binary, 
since the V2 values change significantly by about 150 km s -1 . 
The third emission H a peak is however not clearly visible 
in the GIRAFFE spectra so we can't fully support that hy- 
pothesis by comparing all TODCOR and H a measurements. 
Nevertheless it is worth noticing that during two consecutive 
nights when two R/G observations were taken with roughly 
one hour time span (JD~2454007.5 and 2454008.5) the po- 
sition of the peak on the TODCOR map changed by about 45 
km s _1 . The remaining R/G spectrum (JD~2454006.5) - 
from which the highest value of V2 was derived (116 km s _1 ) 
- had the highest signal to noise ratio. 

From the TODCOR measurements of the B component's 
RVs we derived an acceptable orbital solution with the rms 
of measurements of 7.1 km s _1 and reduced \ 2 °f the fit 
nearly 3. The parameters of this solution are listed in Table 
[3] and the model RV curve, together with the single mea- 
surements and residua (listed in Table [T}, are plotted on the 
right panel of the Fig. [2] In this orbital fit we did not in- 
clude the measurements derived from the H a lines, because 
the resulting parameters were not very different from the 
ones given, and the quality of the fit (in terms of rms and 
X 2 ) was substantially worse. One can see that the systemic 
velocity 7s differs from the 7a value by less than the sum 
of their errorbars, which is expected in case of gravitational 
bounding. This is an argument for the correctness of this 
fit and thus the binary nature of ASAS-01 B, but consider- 
ing the amount and quality of the available data, we are far 
from making any conclusive statements about this hypoth- 
esis. The binarity of ASAS-01 B will be discussed in a later 
part of the paper. 

3.5 H a line profiles 

The existence of strong H a emission is a manifestation of a 
substantial stellar activity in all components of the ASAS-01 
system. The closer inspection of the H a lines reveals their 
probable variability. In Fig. [3] we present H a profiles from all 
of our spectra. On the spectra from HamSpec and UCLES 
all three components are clearly visible and marked. One has 
to remember that for the GIRAFFE setting the H a line was 
close to the edge of the echelle order, so the SNR around 
the line was lower. For every spectrum we measured the H a 
equivalent width and collected them in Table [4] 

The Aa component seems to be particularly strong on 
both GIRAFFE spectra from Sept. 30 (JD~2454008.5). This 
may be due eithe to the real increase of the line intensity 
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Table 4. H a lines equivalent width measurements for all our 
spectra. The measurement errors are estimated to be around 
0.1A. Values for the B component from the GIRAFFE spectra 
(R/G) should be treated with caution. 



JD 


EW Aa 


EW Ab 

[A] 


EW B T/S 


-2450000 


[A] 


[A] 



4006.50332 


0.74 


0.54 


0.18 


R/G 


4007.52761 


1.38 


0.71 


0.15 


R/G 


4007.57128 


1.39 


0.40 


0.13 


R/G 


4008.47697 


1.69 


0.65 


0.10 


R/G 


4008.52042 


2.04 


0.58 


< 0.1 


R/G 


4375.87729 


1.60 


0.67 


0.60 


S/H 


4727.14578 


0.60 


0.25 


0.87 


A/U 



or to blending with the B component. However, in the lat- 
ter case the velocity of B calculated from the line position 
would differ by 50-70 km s _1 from the maximum on the 
coresponding todcor map. In the case of GIRAFFE data 
the intensity and even the existence of the third peak is 
highly dependent on how the spectrum was normalized to 
the continuum. The low SNR of the GIRAFFE spectra and 
the high rotational broadening of the stellar lines make this 
operation quite challenging and its results unreliable, espe- 
cially at the orders' edges (i.e. around the H a line). 

The EW measurements, given in the Table U were 
made under assumption of three emission components for all 
spectra, but we stress that for R/G observation this assump- 
tion may be wrong due to low SNR. This data set however 
shows a systematic decrease of the EW measurements for 
the B component (by a factor of ~ 2), and a rising trend for 
the Aa (by a factor of nearly 3). The EWAb values seem to 
be randomly oscillating around a constant value. Neverthe- 
less these measurements are unsecure and should be treated 
with caution. 

In the HamSpec and UCLES observations the B com- 
ponent is obviously stronger than in those performed with 
GIRAFFE. The Aa component shows values substantially 
different between the two spectra. The EWAb for the Ham- 
Spec spectrum is within the oscillations seen in the GI- 
RAFFE data, but in the UCLES spectrum it reaches the 
smallest value of all data sets. This may however be due to 
the blending with the component B. Given the activity of 
the system, we expect the observed variability of H a to be 
real, however the quality of the data does not allow for a 
firm statement. 

We also want to note that in three of five GIRAFFE 
spectra, as well as in HamSpec and UCLES observations, we 
have found signatures of Hp emission lines. Especially in the 
UCLES spectrum the emission from all three components is 
clear and the EW measurements were possible. The shape of 
the profile is very similar to the H a line. We found EW's to 
be 0.70, 0.17 and 0.62 A for the Aa, Ab and B components 
respectively, with the uncertainty of ~0.1 A. The Hp line 
from the UCLES spectrum is also shown in Fig. [3] 




6560 6565 
Wavelenght 



4860 4865 
Wavelenght 



Figure 3. H a emission profiles variability as seen in our spectra, 
especially for the B component. For the GIRAFFE setup the H a 
line is close to the order's edge, so the SNR around it is very 
poor, which makes the continuum correction and EW calcula- 
tions quite challenging. The Hp line from the UCLES spectrum 
is also shown for comparison. JD-2450000 and spectrographs are 
labelled, components are marked. Spectra were continuum nor- 
malized. 



cated in SAAO, and the 0.41-m Prompt-4 and Prompt- 
5 robotic telescope^, located in the Cerro Tololo Inter- 
American Observatory in Chile. A more detailed descrip- 
tion of the observational settings, reduction procedure and 
calibration to standar d phot omet ric system can be found in 
iHelminiak fc Konackil (|201 lh and iHehniniak et al.l (|2011i > . 

The majority of the SAAO observations were done in 
December 2008. Most of the orbital phases were covered, 
however without the primary minimum in V and with only 
few points in I. A clear out-of-eclipse modulation, probably 
originated by a presence of a cold spot, was noticed. 

The PROMPT observations were carried out mostly in 
July 2009, with two additional nights at the ends of August 
and September the same year. This time several primary 
minima were covered, however more points were collected 
around the secondary minimum. Again we noticed an out- 



4 PHOTOMETRIC OBSERVATIONS 

The photometric observations of ASAS-01 were carried out 
in V and / bands with the 1.0-m Elizabeth telescope, lo- 



3 Panchromatic Robotic Optical Monitoring and Polarime- 
try Telescopes. PROMPT is operated by SKYNET - a dis- 
tributed network of robotic telescopes located around the 
world, dedicated for continues GRB afterglows observations, 
http : //skynet .unc . edu 
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of-eclipse modulation which was bascially the same as for 
the former SAAO data in regards to its shape. 

The last small part of the data, mostly in V, comes 
from SAAO observations conducted in October 2009. Again 
a wide range of the orbital phases was covered, including 
both minima, and an out-of-eclipse modulation was noticed 
to be of the same shape as for previous data sets. This led 
to a conclusion that if the modulation comes from the cold 
spots, the spot pattern did not change much between De- 
cember 2008 and October 2009. 

In total we gathered 598 brightness measurements in V 
and 586 in /, some under quite challenging conditions, which 
gave us a full phase coverage. Several primary and secondary 
minima as well as an out-of-eclipse variation were recorded, 
which allowed us to perform a detailed modeling of the sys- 
tem, starting from an update of the ephemeris and ending 
with spot and third light properties. For the last purpose 
especially useful were SAAO observations from December 
2008, where in good seeing conditions we were able to dis- 
tinguish A and B components on the CCD images (Fig. [T]). 
From a number of images where the components were well 
detached, we estimated their brightness ratio in both bands. 
In this way we deduced that eclipses occur in the brighter 
component A, and obtained reasonable starting values for 
the third-light contribution, which were later incorporated 
in the model. 



5 THE PHYSICAL MODEL 

The modeling procedure was in many aspects identical to 
the ones described in the three previous papers of this se- 
ries. We used the jktebop code (|Southworth et al.l[2004al [bh. 
which i s based on the Eclips i ng Binaries Orbit Program 
(ebop; Popper fc Etzell Il98ll ; lEtzell [l981), and phoebe 
Prsa fc Zwittei 12005]) - a n implementation of the WD code 



Table 5. Absolute physical parameters of the ASAS-01 A eclips- 
ing pair. 



Wilson fc Devinnevlfl97ll ). The jktebop was used on the 



ASAS light curve and served mainly to correct the ASAS 
ephemeris and check for marks of a non-zero eccentricity. We 
did not use it extensively since it does not allow for spots or 
work on RV curves. We should also note that PROMPT 
and SAAO light curves are almost indistinguishable and 
initial fits performed separately gave similar results, thus 
we decided to use all photometric data sets simultaneously 
in the whole fitting process. Photometric data points were 
weighted by their formal errors. Absolute values of several 
physical parameters, like radii, bolomeric magnitudes and 
distance, were obtained with the jktabsdim code, available 
with jktebop, in the same way as decribed in previous pa- 
pers. 

We improved the ephemeris from jktebop with phoebe 
using all available data, and used them to create a final 
orbital fit (Sect. [3]) and later the full solution. We set the 
third light level in both bands to the values obtained from 
the CCD images and let them be fitted for in later stages of 
the process. We consider the uncertainty in the third light 
level as the main source of errors of physical parameters such 
as radii, temperatures and co mponent magnitudes in bot h 
bands (see also ASAS-08 in: iHelminiak fc Konackil l201ll ). 
The third light uncertainty may possibly have its origin in 
the intrinsic photometric variability of component B, which 
is expected for such an active star. 



ASAS-01 A 


Primary 


Secondary 


Parameter 


Value 


± 


Value 


± 


in 




82 ±2 




a [R ] 




2.502 ± 0.035 




M [M ] 


0.612 


0.030 


0.445 


0.019 


n 


4.93 


0.13 


5.21 


0.19 


R [R Q ] 


0.596 


0.020 


0.445 


0.024 




4.68 


0.31 


4.79 


0.45 


Vrot [km s _1 ] a 


67.3 


2.5 


50.5 


2.8 


V -I [mag] 6 


1.84 


0.14 


2.17 


0.14 


Teff [K] 


3750 


250 


3085 


300 


M boi [mag] 


7.76 


0.30 


9.23 


0.45 


M v [mag] 


9.18 


0.31 


10.65 


0.89 


d [pc] 




39 ± 


6 




Parameters 


of spots 






Longitude [rad] c 


0.20 


0.05 


4.71 


0.07 


Radius [rad] 


0.7 


0.1 


0.3 


0.08 


Contrast 


0.9 


0.1 


0.9 


0.1 



a In case of a synchronous rotation 

6 Colour dereddened with E(V - I) = 0.02 mag 

c Counted counter-clockwise with defined as the direction 

toward the companion star 

The starting values of effective temperatures were de- 
duced from the empir ical colour-temperature relations by 
IWorthev fc Leel |201ll ) on the basis of dereddened V — I 
values for every component. For the dereddening we used 
E(B-V) = 0.0145 value from lSchlegel et all l| 19981 ). In gen- 
eral phoebe allows for calculating fluxes in every passband 
separately for every star. If one has at least two light curves 
in various bands, the program also allows for fitting both ef- 
fective temperatures simultaneo usly using the concep t of a 
"binary effective temperature" l|Prsa fc Zwitterl 120051). The 
treatment of limb darkening and reflection albedos was the 
same as in our previous papers, as well as the way of finding 
parameters of cold spots. We tested several different config- 
urations and found the best solution is when a large cold 
spot is located on the secondary and is obscured during the 
secondary eclipse, and a second small spot with maximum 
appearance around phase <f> = 0.25 on the primary. Because 
the light curves do not contain much information about the 
latitude of spots, we decided to locate spots near the stars' 
equators and keep this location fixed. Moving towards other 
latitudes did not change the final model significantly - the 
difference in inferred values of stellar parameters was smaller 
than uncertainties coming from other sources. 

The final absolute values of derived physical parameters 
of ASAS-01 A are collected in Tabled In Fig. 0] we present 
the observational and modeled light curves in V (left) and 
I (right) bands. One can see that the data taken with two 
different instruments and within almost a year overlap quite 
well, so it was justified to use them simultaneously. As was 
mentioned before, the pattern of the spots did not change 
significantly during this time, and no short-scale evolution 
occurred, as was found for two other low-ma ss systems - 
ASAS-09 and ASAS-21 (IHelminiak et alj|201ll ). The 3D re- 
construction of the eclipsing pair in three orbital phases 
- 0.04 (primary eclipse), 0.25 (quadrature), and 0.46 (sec- 
ondary eclipse) - is presented in Fig. [5] (top and middle), 
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Phase 



Figure 4. The observed V (left) and / band (right) light curves of ASAS-01 obtained from our SAAO (+) and PROMPT ( X ) observations. 
An overplotted red line shows the best-fitting model. An out-of-eclipse light modulation, coming from the presence of two cold spots, 
is clearly visible. A significant change of the difference between the depths of the minima in V and / bands indicate the temperature 
ratio is substantially different from unity. Residua are plotted on the lower panels. Colour figure is available in the on-line version of the 
manuscript. 



together with the influence of spots on the light curve (bot- 
tom left) and radial velocities (bottom right). On the ARV 
plot we also present RV residuals in the same form as in Fig. 
[2] One can see that the influence of spots on the RV curves 
is much smaller than the spread of measurements, which al- 
lows us to conclude that the uncertainty of the mass is not 
underestimated. We did not succeed in reaching a 3% level 
of precision in mass and radius determination but we were 
close (3.4 - 5.4%). 

One can also note that our distance determination - 
39(6) pc - is in very good agreement with the result ob- 
tained by iRiaz et aL (120061). Assuming the angular separa- 
tion of 1.405(3) asec ( Bergfors et aljfeoicj ). the projected 
physical separation of the A and B components is thus 
a.AB = 55(8) AU. This corresponds to an orbital period 
Pab — 333 yr, assuming a circular orbit and the total mass 
of component B Mb — 0.45 Mq (similar to Ma6 because 
of almost exactly the same value of V — I), or 292 yr for 
Mb — 0.9 Mq (see Sect. 16.21 for a discussion). The gravita- 
tional influence of B on the pair A would not be detected 
with eclipse timing or precise RV measurements; however 
the orbital motion of B should be detectable with adaptive 
optics or various interferometric techniques. 



6 DISCUSION 

6.1 Evolutionary status 

We compare our results with three se t s of theoretical 
isoch rones: Yonsei-Yal e (Y 2 ; Yi et al. 2001; Demarque et al.l 
l2004h, Dartmouth dDotter et all l2007t) . and Padova 
( Girardi et "aT]|2000l ; [Marigo et al.ll2008l l. 



The comparison is shown in Figure [6] We present the 
bolometric magnitude, radius, effective temperature, abso- 
lute V magnitude and the V — I colour as a function of stellar 
mass. We compare our results with isochrones for an age of 
1 Gyr and two cases of metallicities: (1) solar, Z ~ 0.02 (red 
lines), and (2) above solar (green lines), however different 
for every set: Z = 0.04 for Y 2 , 0.035 for Dartmouth, and 
0.03 for Padova. Since the accuracy of our measurements is 
worse than the canonical 3%, the isochrones are shown only 
for comparison, i.e. we do not attempt to determine the age 
of the system. It is however enough to conclude that ASAS- 
01 is probably a main sequence object. The metallicity de- 
termination is even more insecure, however from the Fig. [6] 
one can claim that metallicities higher than solar are pref- 
ered - on the M vs. Mboi and logT e // planes the Z > 0.02 
isochrones fit significantly better to the data points. 

The main sequence evolutionary stage of ASAS-01 is 
confirmed by its galactic kinematics. We used our deter- 
minations of the systemic velocity and distance together 
with the position and proper motion from the PPMXL 
catalogue: ft a _ 120.4 ± 4.0 m as yr -1 , jj,s = —36.7 ± 
4.0 mas yr" _1 (|Roeser et al.ll20ld) . The obtained values of 
U = 3.3 ± 2.3 km s _1 , V = -16.8 ± 2.1 km s" 1 and 
W = -7.8 ± 2.3 km s _1 put ASAS-01 in the galactic thin 
disk, and only margi nally in one of the young moving groups 
recently reported bv lZhao et al. (2009, ID=15 therein). 

It is worth noticing that on the M/R plane the measure- 
ments agree within errors with the Y 2 and Dartmouth mod- 
els. The second one especially seems to reproduce the physi- 
cal properties of ASAS-01 A components. It is in a somewhat 
contrary to many recent r esults of low-mass d etached eclips- 
ing binaries studies (see: iKraus et al.l l20l"ll ). The usually, 
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Figure 5. Top and middle: Three-dimensional reconstruction of the ASAS-01 A eclipsing binary, based on our best-fitting model in 
quadrature (top) and during eclipses (middle). Dimentions are scaled so that a = 1. Bottom: Influence of spots on the light curves (left) 
in V (solid line) and / (dashed) bands, and on the RV curves of the primary (solid) and secondary (dashed) component. On the RV 
panel we also plot the residua of the model, in the same way as in Fig. [2] On the LC panel, the large scatter on the V curve around 
(p ~ comes from the sampling of the stellar surface used in the PHOEBE code. 



but not always, observed characteristic of low-mass stars in 
close binaries is that theoretically predicted radii are smaller 
and temperatures larger with respect to what is observed 
(Ribas et al. 200||), an d this situation is seen for the Padova 
models. The Dartmouth models also seem to follow that 
trend, but only by a few per cent. Within our uncertainties 
they correctly predict the observed radii, and, in the case of 



Z — 0.035, also temperatures. The deviation is smaller for 
the larger component, which presumably rotates faster, thus 
it does not support the possible radiu s-rotation relation for 
M < 0.7 M Q stars l|Kraus et al.ll201ll ). At the same time the 
properties of the secondary are very well reproduced by the 
Y 2 models (except for the colour) but for the primary we 
see the "typical" underestimation of the radius and overes- 
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Figure 6. Comparison of our results for ASAS-01 A eclipsing 
pair with theoretical 1 gyr isochrones on mass M vs. (from top 
to bottom): bolometric magnitude, radius, effective temperature, 
absolute V magnitude, and V — I colour. Y 2 isochrones are de- 
picted with solid, Darthmouth with dashed, and Padova with 
dot-dashed line. Isochrones for solar metallicity (Z ~ 0.02) are 
depicted with red and for higher Z with green. The colour ver- 
sion of the picture is available in the on-line version of the paper. 



Table 6. Photometric properties of the component B from the 
PHOEBE model 



Parameter 


Value 


± 


Fv,3 W°\ a 


34.8 


1.5 


Fx, 3 [%r 


42.9 


3.2 


(V - 7) 3 [mag] 6 


2.15 


0.13 


M V}3 [mag] 


10.82 


0.70 



a Fractional fluxes defined as a percentage of the flux of 
component A 

b Colour dereddened with E(V - I) = 0.02 mag 

Table 7. Values of the inclination of the component B putative 
orbit for various values of qg and masses of the more massive 
component of 0.45 Mq (ig 0.45) an d 0.5 Mq (ig o.s)- The mass 
function is taken from the Table [3] 





M Ba sin 3 ig 
[Mq] 


M Bb sin 3 ig 

[Mq] 


[°] 


»S,0.5 
[°] 


1.00 


0.226 


0.226 


53 


50 


0.95 


0.251 


0.238 


55 


53 


0.90 


0.280 


0.252 


59 


56 


0.85 


0.315 


0.268 


63 


59 


0.80 


0.358 


0.286 


68 


63 


0.718 


0.45 


0.323 


90 


75 


0.684 


0.5 


0.342 




90 



/(Mb.) = 0.0565(92) Mq 



timation of the temperature. From this short discussion we 
can only conclude that the Dartmouth models may be the 
best to reproduce the properties of low-mass stars (see also: 
iThompson et al]|2010h although discrepancies at the level of 
~3% are present. This is however significantly smaller than 
for the majority of previous studies where 5-15% difference 
in radii were claimed. 

6.2 The companion B 

Having the two components separated on V and / band im- 
ages (Fig. [TJ, and the physical model created with phoebe 
with the third light included, we were able to estimate the 
photometric properties of component B. We present the frac- 
tional fluxes in both bands, the dereddened V — I colour, 
and the absolute V magnitude in Table [6] 

In Figure [7] we present our measurements on the 
V — 1 1 My plane. Data from Tables [5] and [6] and the same 
isochrones as in Fig. [6] are used. We see that again the Dart- 
mouth set for Z = 0.02 gives the best match for the eclipsing 
pair A. Component B is about 0.7 mag brighter but has al- 
most the same colour as Ab (V — I — 2.14 mag). Due to 
relatively large errorbars (magnified mainly by the uncer- 
tainty in the third light parameter in our model) one can 
formally find V magnitudes of B and Ab almost equal and 
still rather consistent with the considered Dartmouth and 
both Y 2 isochrones. However, if assuming that component 
B is composed of two nearly-twin stars with masses ~0.45- 
0.5 Mq, which would thus have the V — I colour around 
2.15 mag, i.e. almost the same as the value observed for B, 
those putative Ba and Bb components would be located in 
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Figure 7. Colour-brightness diagram with our measurements for 
ASAS-01 Aa, Ab and B components, together with theoretical 
1 Gyr isochrones. Line and colour coding is the same as in Fig. [6] 
The square shows the approximate position of the components of 
ASAS-01 B if it was composed of two identical stars showing the 
same V — I value as observed. The colour version of the Figure is 
available in the on-line version of the manuscript. 



Fig. [7] on the position marked by the empty square (low- 
ering the flux by a factor of 2 increases the magnitude by 
~0.7). In such a situation, the Darthmouth Z = 0.02 and 
Y 2 isochrones would reproduce the observed properties of 
the system much better. We consider this fact as a support 
for the hypothesis of the binarity of ASAS-01 B. 

The binarity of component B, with a mass ratio close 
to 1 and component colours close to the value observed, is 
consistent with the fact that no marks of additional eclipses 
are found in our V and I band light curves. From the orbital 
fit for component B radial velocities, presented in Table[3] we 
can estimate the orbital inclination is for various values of 
the mass-ratio qs and the desired mass of the more massive 
component Msa- We present those calculations in Table [7] 

For a number of different values of qs, we present val- 
ues of M sin 3 ib for Ba and Bb, and values of inclination 
angles, for which Mbo. = 0.45 or 0.5 M . From this Table 
one can see, that if ASAS-01 B is a binary, to explain all the 
observed properties one needs the inclination angle is to be 
between 50 and 70 degrees, which corresponds to q ^ 0.8. 
For these values of masses of Ba and Bb we can expect the 
observed colour to be around the value given in Tab. [6] As- 
suming lower values of qB we would end up in a situation 
where additional eclipses occur and the Bb component is 
significantly fainter. In such a case the total brightness of 
ASAS-01 B would be lower than it is observed, however still 
within relatively large error bars. 

If ASAS-01 is a double-double system, it would join a 
small group of this kind of interesting objects. To date only 
two other low-mass detached eclipsing binaries are known 
to be in similar configuration: BD -22 5866 - a system 
with K7+K7 eclipsing pair and M1+M2 non-eclipsing bi- 
nary (|Shkolnik et al.1 l2008h . and YY Gem - the faintest 
member of a sextuple system composed of three spec- 
troscopic binaries - Castor A, B and C (a Gem ABC ; 
IVinter Hansen. Neubauer fc Roosen-Raadlll94cl ; iKronll 19521 ; 



iBopd 119741 ). There are also few examples of LMDEBs 
known to have a sin gle ad ditional companio n, like LP 133- 
373 (|Vaccaro et alj l2007h. HIP 965lfl (iHuelamo et af] 
l2009h . MR Del JPribulla et alJl2009l ; lDiurasevic et aljfeoilf). 
NLT T 41135 Jlrwin et al l l201Ch ASAS-08 (|Montes et all 
120071 ; iHelminiak fc Konackil l201ll ') or the triply-eclipsing 
KOI- 126 ( Carter et al.ll201lr ). Such systems not only allow 
more rigid constrains on the evolutionary models than in 
the cases of "lonely" eclipsing binaries, but also play im- 
portant role in testing star-formation theories, stellar pop- 
ulation codes and dynamical interactions in multiple stel- 
lar systems. The relative brightness and small distance to 
ASAS-01 make it a valuable object for further studies. 

7 SUMMARY 

As a result of our analysis of the components of 
ASAS J01328-3821.1, an M-type multiple system, we ob- 
tained a complete set of orbital and physical parameters of 
the detached eclipsing pair A, and photometric and spec- 
troscopic properties of component B. Due to a number of 
difficulties our results are of an accuracy of 3-6 %, which is 
not enough for performing reliable estimation of the evolu- 
tionary status of the system. Within the errorbars we may 
however reproduce our results with main sequence theoret- 
ical models. The observed spectroscopic and photometric 
properties of component B of the system suggest that this 
is a binary composed of nearly identical stars. However the 
eventual orbital solution is quite uncertain, and Mv and 
V — I errors are quite large, thus we can not make any final 
conclusions on this topic. ASAS-01, as a member of a rare 
class of late-type multiples with eclipsing detached compo- 
nents, should be considered as a potential testbed for stel- 
lar formation and evolution theories and definitely deserves 
more attention and investigation, especially considering the 
system's spectroscopy and photometric properties of com- 
ponent B. 
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